Coherent response of a low T c Josephson junction to an ultrafast laser pulse 
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By irradiating with a single ultrafast laser pulse a superconducting electrode of a Josephson 
junction it is possible to drive the quasiparticles (qp's) distribution strongly out of equilibrium. 
The behavior of the Josephson device can, thus, be modified on a fast time scale, shorter than the 
qp's relaxation time. This could be very useful, in that it allows fast control of Josephson charge 
qubits and, in general, of all Josephson devices. If the energy released to the top layer contact SI 
of the junction is of the order of ~ /ij, the coherence is not degradated, because the perturbation 
is very fast. Within the framework of the quasiclassical Keldysh Green's function theory, we find 
that the order parameter of SI decreases. We study the perturbed dynamics of the junction, when 
the current bias is close to the critical current, by integrating numerically its classical equation of 
motion. The optical ultrafast pulse can produce switchings of the junction from the Josephson state 
to the voltage state. The switches can be controlled by tuning the laser light intensity and the pulse 
duration of the Josephson junction. 

PACS numbers: 74.50.+r, 74.40.+k,74.25.Gz, 74.25. Fy 



I. INTRODUCTION 



The characteristic frequency in the dynamics of a Josephson junction (JJ) is the so-called Josephson plasma fre- 
quency ujpj (e.g. 10 -j- 100GHz). Coupling of a JJ to a microwave field leads to the well known lock-in conditions, 
which show up as Shapiro steps in the I/V characteristic. On the other hand, photo-response to radiation in a su- 
perconductor induces heat relaxation (bolometric effect) and non equilibrium generation of quasiparticles (qp) 2 *^. 
Both phenomena are extensively studied since they are relevant for the fabrication of fast and sensitive detectors. The 
models used are phenomenological^^*^, mainly involving different temperatures associated to separate distributions 



of electrons and phonons out of equilibrium. 
Recently, laser light with pulses of femtosecond duration r c 



€ (10 14 s, 10 13 s) has become available, as a source to 
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FIG. 1: Sketch of the Josephson junction exposed to a laser radiation pulse. 



test the photo-response of a Jj£. Ultrafast pulses can be extremely useful, in that they allow studying an unexplored 
regime in non equilibrium superconductivity. Indeed, photon absorption, by creating electron-hole (e-h) pairs at very 
high energies, drives the quasi-particle (qp) energy distribution out of equilibrium during the time r c . The qp non- 
equilibrium distribution depends on the energy relaxation time parameter te, defined as the time by which a 'hot' 
electron is thermalized by repeated scatterings with other electrons or phonons. The process involves generation of 
many qp's during energy degradation, until the system relaxes back to the equilibrium distribution function n (u). 



This time scale is determined by the electron-electron interaction time r e _ e and the electron-phonon interaction time 
T e -ph, which are strongly material dependent 9 , ranging from 4 • 10~ 7 s for Al to 1.510 _10 s for Nb. In this work we 
analyze the possibility that, keeping temperature quite low, ultrashort laser radiation induces direct switches out of 
the Josephson conduction state at zero voltage, due to coherent reduction of the critical current J c . 

There are many reasons for the switching from the zero to the resistive state in a Josephson junction. Among 
these, thermal escaped, quantum tunneling 11 , latching logic circuits^ and pulsed assisted escaped. A clear cut 
discrimination between different mechanisms can be difficult to achieve. In our case quantum escape is ruled out 
because the temperature is not expected to be low enough. Also, we assume that there is no external circuit to induce 
switching and re-set of the zero voltage state as in latching logic elements. 

Pulsed assisted escape is a generic term for a large class of phenomena including in principle bolometric heating of 
the junction which is re-set in relatively slow tiniest. Production of quasiparticles generated by X-ray radiation has 
been studied up to recentlyi^i^. A cascade follows, which increases the number of excitations and lowers their energy 
down to the typical phonon energy lojj in a duration time, which is of the order of the nanoseconds. Subsequently, 
qp's decay by heating the sample. However, the power of the laser can be reduced enough and both the substrate 
and the geometry can be chosen such that the energy released by the radiation on the junction can be small. On the 
other hand, appropriate experimental conditions can make the time interval between two pulses long enough, so that 
the bolometric response is negligible. 

Generally speaking, junctions are more sensitive to pulses especially when their harmonic content is close to lo p j, 
but this is not our case. In fact the laser carrier frequency (Q ~ lOOTHz) is quite high compared to u> p j and we 
consider the case Q >> t~ x > ui p j > t^ 1 , what implies that little relaxation takes place during the duration of 
the pulse. t c should also be shorter than the pair-breaking time h/A Q ~ 1 bps. Here A Q is the unperturbed gap 
parameter. Our approach assumes that, on a time scale intermediate between the pulse duration and the relaxation 
time ~ te, the order parameter of the irradiated superconductor is sensitive to the non-equilibrium qp distribution, 
which modulates it coherently till it switches out of the zero voltage state. 

To analyze the dynamics of the order parameter and the way how the latter affects the Josephson current, we 
adopt a non-equilibrium formalism based on quasiclassical Green's functionsiSiii. The quasiclassical approach has 
been mostly used in the past in connection with the proximity effect^, as well as with non-equilibrium due to other 
space inhomogeneity conditions^. As far as we know, this is the first time that its extension to non-equilibrium in 
time is applied to a coherent response after an ultrafast laser irradiation. 

The quasiclassical approximation to the Gorkov equations, is obtained by averaging over the period of the optical 
frequency f2, which is a fast time scale 19 . Our equations include the physics of the cascade process, which occurs when 
one focuses on the kinetics of the qp diffusion. A kinetic equation approach to the steady state non equilibrium qp 
distribution, including phonon scattering has been developed in ref^*, for light irradiation, mostly in the microwave 
range. The cascade regime is extensively discussed in Ref. 15 , however it will not be specifically addressed here. 

Instead, if the switching of the irradiated superconductor due to the ultrafast pulse takes place prior to the occurrence 
of qp relaxation, an approximate solution of the dynamical equations can be derived, which describes an istantaneous 
response of the order parameter. 

We take a low T c JJ with an s-wave order parameter as the reference case (e.g. a high quality Nb or Al junction) 
and T << T c . The optical penetration depth of the laser light \$ in the topmost superconductor exposed to radiation 
SI is assumed to be shorter than its thickness, so that any modification induced by the radiation field only involves 
SI itself (see. Fig.JU). In a small size J J the spatial variation of the order parameter along the lateral dimension of 
51 is not taken into account, except when the qp diffusion process cannot be ignored. 

We consider just one pulse of given duration t c which releases the energy £ per pulse, by exciting e — h pairs 
and by creating a non equilibrium distribution of qp's. A related dimensionless quantity q, as defined in Eq. 10, 
parametrizes the strength of the perturbation due to the radiation. The perturbation is assumed to be small so that 
only the lowest order in the expansion in q is retained. This allows us to derive a temporary reduction of the order 
parameter A induced by the pulse, as shown in Fig.©. We do not give a detailed description for the relaxation of the 
non-thermal qp distribution in the irradiated superconductor. The self-energy terms corresponding to this process 
require further analysis. According to the Eliashberg formulation 9 , these terms affect the quasiparticle amplitude Z(u>) 
introducing changes in the phonon distribution and retardation in the response. Nevertheless, we expect that these 
self-energy terms become effective only on a longer time scale after the laser pulse. Our equations pinpoint a non 
retarded evolution of the order parameter prior to relaxation, which implies a reduction of the critical current. This 
shows that the coherent modulation of the gap parameter can produce switching of the junction out of the Josephson 
state. 

The switches are studied numerically by solving the classical equation fo motion of a current biased JJ with current 
J close to the critical curent J c , during the excitation process. After the switching the dissipation in not treated 
selfconsistently: a standard dissipation, typical of thermal equilibrium, is assumed in the JJ dynamics, by adding 
a conductance term in the numerical simulation. We stress that the assumed model for dissipation determines the 



actual qp branch of the I-V characteristics, but does not affect substantially the switching probability. The switching 
from the Josephson state to the voltage state in the parameter space (q,r c , J/J c ) is reported in FigEl Interestingly, 
we find that for fixed value of q and J / J c there is an optimum pulse duration to achieve the strongest sensitivity of 
the junction to the switching process. We show that this is due to the way how the non equilibrium qp's distribution 
affects the pairing in SI. The paper is organized as follows. 

In Section 2 we calculate the non equilibrium qp distributions immediately after the pulse, prior to relaxation. 

In Section 3 we introduce the time dependent quasiclassical Keldysh Green functions formalism extended to the time 
domain, using the general frame given in appendix A. We calculate the correction to the single particle propagators 
to first order in q. A correction to the gap and to the Josephson critical current follows due to the laser pulse. 

In Section 4 the dynamics of the JJ after the pulse is simulated numerically. Finally, a summary of our results is 
given in Section 5. 

Appendix A collects the formulae of the quasiclassical approximation in the non equilibrium Keldysh theory, which 
are used in the core of the paper. In Appendix B we derive the kinetic equation for the non equilibrium distribution 
function which drives the relaxation of the system. The equations of motion for the quasiclassical retarded Green 
function is reported in appendix C, while the equation of motion for the advanced and Keldysh Green function can 
be derived in the same way . 

II. THE NON-EQUILIBRIUM QP DISTRIBUTION 
A. Non-equilibrium electron-hole pair excitations induced by optical irradiation 

The optical frequencies (0 ~ lOOTHz) building the wavepacket of the laser pulse excite e — h pairs at high 
energies. As explained in the introduction the non-equilibrium arising from the alteration of the qp distribution has a 
relaxation time te which is long compared to the optical period: Qte >> 1. In addition to this, the duration of the 
pulse r c ~ w" 1 is even shorter than the pair breaking time, so that we expect that, in our case, dissipative phenomena 
do not affect the coherence of the superconductor on the time scale r c . 

Qp's are generated as if the metal were normal, because superconductivity doesn't play any role in their excitation 
at large energies. They propagate according to something very much like the free particle time-ordered Green's 
function g Q (from now on we put % = 1): 

g T {k- t-t')= -i(T[c k (t)cl(t')]) = -ie~^ 1 ^ [(1 - n k )0(t - if) - n k 9(-t + if)] . (1) 

Here £ is the qp energy with momentum k and is measured from the chemical potential fi. n k is the qp distribution 
function. We assume that e—h symmetry is conserved in the excitation process, so that fi is not altered with respect 
to its equilibrium value. 

The equation of motion for the Green's function g in the presence of the radiation field is: 

(i^ t -^[V r -^A(r,R,t)] 2 + v\ ~g(?,R,t,t') = 5(r)8(t-t'), (2) 

f is the relative space coordinate, while R is the center of mass coordinate. The vector potential is a wave-packet 
centered at frequency f2 according to: 

A(r, R,t)=^J2^ ® £ t)e^-nt) . (3) 
± p 

Here a± are slowly varying 'envelope' functions of R on the size of the irradiated spot and on the time scale SI -1 . We 
look for solutions of eq. (J2J in the form: 

g( f, R; t, t') = g{r, R; t, t') + £ £ g ± (p, R, t, <>^"-o t ) ^ (4) 

± p 

where g and g ± are slowly varying functions of R and t on the same scales. A similar expansion can be done w.r.to the 
variable t' . Following Eq. J3J|, a decomposition of eq.(j2J) into harmonics arises^. We define the zero order harmonic 
equation as the one that does not contain exponentials e ±lQt . By averaging over a period SI -1 we neglect harmonics 



of order two, or higher. This amounts to include one photon excitation processes only, with released energy E . Some 
extra details can be found in Appendix A: 

d 1 



dt 2m mc- 



" v ' + — E S + (0,R,t)a^',R,t)e^-p'y p + M g(r,R,t,t') = 5(r)S(t - t') . 

f ,p" / 

Fourier transforming w.r.to f ( r — > p — > £) we have: 

g{p, R, t, + -^E S + % (p"' ^ f ) - P" - # # *. = *(* ~ ■ ( 5 ) 



mr 



The radiation field generates and annihilates high energy e — h pairs. Hence we assume that the forcing term conserves 
the total impulse, p 1 + p" ~ 0, but p 1 — p" transfers an energy 2£ to the electrons. Therefore we take the coupling 
term in the Hamiltonian as: 



'7T 



A Gaussian shaped time dependence has been chosen for the pulse with half- width u> c , while the space dependence 
has been neglected for simplicity. In Eq. © the dimensionless quantity appears: 



q ~ c n mfi 2 i?2 ' 1 j 

Where R is the laser spot (see Ea 11611 . Here the number of excited e — h pairs is ~ SI/u>d, with uj d the Debye energy. 
Experiments 2 *^ show that the energy released by the pulse can be very low, so that we will always expand in q. In 
fact, while in the case of an rf radiation q ~ 1, in the case of a femtosecond laser pulse q ~ 0.01 -j- 0.1, corresponding 
to a fraction of [iJ released per pulse on the superconducting surface of ~ 100/jto 2 . 
The zero order harmonic equation, Eq. JSJ, becomes : 

(idt - <?(£; t, + q^e~^ t2 g(-t; t, f) = S(t - t')) . (8) 

\/7T 



To derive the non-equilibrium correction to the qp distribution function, the kinetic equation IB4I should be solved. 
In place of this we proceed in this work in an heuristic way. Our approach lacks mathematical rigour, but singles out 
directly the role of the laser induced e — h excitations at frequencies (ft — oj c , Q + u> c ). Our result is valid in the limit 
of large £'s and zero temperature, before relaxation takes place. 

We solve Eq© for the retarded Green's function for t > and —t' ~ + by truncating the Dyson equation to 
lowest order in q: 

9 R & t, f ) = t-t')- q ^=J dt"g*& t, fV^sfH; *" - 

= g R (C;t-t')+i q ^% / p 6 , . n+ / dte^K-i^ , (9) 

where g R (ui) = {oj — £ + i0 + }~ 1 is the Fourier transform of the retarded Green's function. The time integral can be 
expressed in terms of the function w[z] = e~ z erfc (—iz). If we now approximate Eq. by evaluating w[z] only at 
the pole and use the the integral representation of the step function: 



9(t) = — / du for Quiz < , (10) 

2m J_ QO uj - z 

the correction 5g R (^,t) to g R for £ >> 0, which includes the non equilibrium qp 's distribution, is: 

8g R {^t,Q-) = -^=e-^- i0+ > «, [-2£/(V2w c )] 9{t), (11) 
From ea. (|ll|l we obtain the time ordered Green's function for t > 0, t' = 0~ and £ > 0: 

g& t > 0, 0") = (_i) e -<«-«> + >* 9{t) (l - -L p [ji£/(V2Wc)] ) , (12) 



with 



p[a 



ds i 



(13) 



p[x] increases linearly with x and it decreases slowly, as 1/x, at large arguments. In Eq. i|12|) we have neglected 9?e[iw] 
because |f| >> 0. 

Eq. H12|l is to be compared with the free propagating time ordered Green's function of eq.Q for the same time 
arguments. Comparison yields the amount by which the distribution function is driven out of the equilibrium: 



Sn(0 » -j= p 



2Z/{V2lu c )] for |f | » A 



(14) 



Note that the expression of Eq. (|14|) changes sign according to sign {£ } . This stems from the assumed e—h symmetry. 
In turn this implies that no charge imbalance occurs. 

Eq. I|14|l can be considered as the non equilibrium distribution for qp 's starting at the time of the pulse t ~ + . 

In the absence of relaxation, a change in the available qp density of states follows. Because (g R (£,,oj))* = g A (t;,uj) 
if u) is real, the correction to the density of state Sv{^ ) is: 



5v(0 = — (6g R - 8g A ) (£,0+) 



2£/(V2w c ) 



for f > 



(15) 



The first stages of the relaxation process involve the inelastic diffusion of qp's in the medium which is qualitatively 
discussed in the next section. 



B. Inelastic diffusion of the qp's at initial times 

Let us discuss shortly what was neglected in the derivation of the change in the equilibrium qp distribution Sn(£) 
given by Ea. lCHfl . 

The single particle Green's function g(p, R, t, £') is assumed to be a slowly varying function of (t + t')/2 = t and a fast 
varying function of t — if . Fourier transforming w.r.to the latter variable ( see Appendix A ) there is an w dependence 
even in the stationary case ( i.e. with no t dependence). This u> dependence is determined by the frequency dependent 
Eliashberg e — ph coupling a 2 F(uj) and is contained in the e — ph self-energy M e - p h(u})2Q. Accordingly, the complex 
qp renormalization parameter Z n (u>) is defined by [1 — Z n (uj)]u) — M e - p h(uj). In our derivation we have not included 
the self-energy, so that we are implicitly taking Z n (u>) — > 1, what applies for large a>(~ f2), prior to relaxation. 

Moreover, because the system is in the superconducting state, we should have dealt with the corresponding 
superconducting parameter Z s (u>). The latter is derived together with the complex gap parameter A(u>,t) with 
A(A a ,t = 0~) = A D from the coupled Eliashberg equations ( we drop the overline on t in the following ). 

The procedure of averaging over the fast time scale fl^ 1 singles out two frequency components of A(cj, t) and Z s (u>, t): 
lu = A a and u> = £1 as a consequence retardation arises from frequencies up to IOlud is neglected, is so large that 
Z 8 and Z n do not differ sizeably. In fact, their real parts differ by a quantity of the order of (ujdA /VL 2 ) 2 \jx{ojd/ A q ). 
A(f2, t) itself is expected to be so small that it can be neglected alltogether. Indeed, in connection with eq. IA10I of 
Appendix A we do not discuss the self-energy terms. Of course this approximation breaks down on the time scale of 
e — ph relaxation. 

Let us now discuss the t dependence. The equation of motion for the qp distribution function n(i?, t) is derived in 
Appendix B, where we take hit — 0, because we neglect charge imbalance corrections. 

In averaging over a few optical periods the kinetic equation for Shl, the electric field E averages to zero. The 
qp relaxation is governed by the collision integral I[n(R,t);t] which describes the inelastic processes. In Ref. 1 — the 
cascade of the e — h excitations due to inelastic scattering is studied in detail. Two stages occur. In the first stage 
e-e interactions multiply the number of excited qp's in the energy range from Q to which is taken as the cutoff 
energy of the pairing interaction. This happens in a time interval short w.r.to the pulse duration (~ 10~ 14 s). In the 
second stage a much slower relaxation process takes place, by which the energy of the qp's reaches A. This process 
involves electron-phonon scattering on a time scale hcu^/A 3 ~ 10 2 ns which is much larger than any time scale in our 
problem. Here we will leave this stage aside. In the time interval we are concerned with, we have little relaxation and 
the energies involved are to » A. 

According to Eq. (|B2|I the distribution function prior to relaxation deviates from the equilibrium value by the 
quantity Sn^ = — 26n(£) given by Eq. (jl4(l . There is no explicit dependence on u> in our correction, becauise 



retardation is nelgected. Still qp's diffuse in space inside the junction over a characteristic distance R ~ y/ 'Dr e - e , 
where D is the diffusion coefficient. Hence 



Sn L (tR,t) = -2 1 ^ Wt Sn(0e{ ^¥1 . (16) 

For relatively large times T e - p h >> t » R^/D we will ignore the spatial dependence by putting R = 0. This is the 
first step of a perturbative analysis of the non equilibrium distribution functions. 

III. CHANGES OF THE SUPERCONDUCTIVE PROPERTIES ON THE TIME SCALE u' 1 

A. The correction to the gap parameter 

In this Subsection we derive the Keldysh Green's function in the presence both of a time dependent gap A(i) and of 
a non-equilibrium qp distribution as given by Eq. I|16() . We assume weak coupling superconductivity and we neglect 
here the frequency dependence of the e — ph coupling parameter A. This follows from the neglecting the retardation 
effects mentioned in Sec. lib on time scales much faster than the e — ph relaxation time. From the Keldysh Green's 
function g K (where the hat denotes matrix representation in the Nambu space, see appendix A) we recalculate the 
gap self-consistently, according to the formula: 

m = [°° < f*ffl\p\, u ,t) > V ~ F du; . (17) 

J — oo 

The average over the direction of the momenta on the Fermi surface is indicated. The Keldysh Green's function in 
thermal equilibrium is: 

gf=tanh^(^-^) . (18) 

Out of equilibrium we use the definition: 

g K = g R h-hg A . (19) 

However, h defined in appendix B is here diagonal, because we assume that no charge imbalance arises. Hence, up to 
first order in q, 

Sg K « n°M - g£ d ) + 6n L (g« - g A ) . (20) 

Here n° L = tanh(/3w/2) is the equilibrium distribution and g^ d — g A d is introduced in Appendix C (see Ea lG2|l and is 
discussed in the following. 

We now first derive the contribution coming from the second term of Eq. 112011 . We start from the outset using Eq. 
(|16[) and performing the quasiclassical approximation. The latter involves an energy integration: 

■ <.+oo 

Sn L • {g« -9f) = - d£ 6n L (Z, t) ■ w, t) - g A (^ u, *)) . (21) 

^ J — OG 

Using the equilibrium BCS functional forms, the Green's functions appearing on the diagonal of g A ' R areSi: 

7 ? 
U V 

^'"^ v-bU^ * + % + (22) 

g A (t w) = * . n+ + _i_ p - n + ' ( 23 ) 



The equilibrium values for and are: 



ki+bY , «g = fi(i-i)] . (24, 



2' E' ' * \2 y E 



with E = a/£ 2 + |A G | 2 . From now on we will drop the subscript in the equilibrium gap parameter (i.e. A = A Q if no 
time dependence is indicated explicitely) . Because the factor <5nL(£,i) appearing in ea. (|21|l . as given by Eq. 1)16(1 is 
odd w.r. to £ only the second term in u 2 and v 2 survives, when the integral in Eq. I(21|l is performed. Let us consider 
the case to > A only and specialize Eq. (|21[) to its diagonal part. According to Eq. I|16[) we have: 



9£e {Sn L ■ g*} = 9fc 



y/2i 



+ OC 



df— p 
^2E H 



uj-E + i,0+ uj + E + iO+ 
= -q(t) 2V2p \2(uj 2 - \A\ 2 )^/{V2uj c ) 



Here we have defined the function q(t): 



(25) 



q(t) = q 



nRi 



Rl + Dt 



(26) 



Doing similarly for g A and subtracting, the imaginary part cancels: 

Sn L • (g? ~ 9*) = -q(t)4V2 p \2{oj 2 - |A| 2 ) * /{y/2u c ) 



, for to > A 



(27) 



Here the largest contribution of the non-equilibrium excitations arises from to ~ uj c . On the other hand lo c can be 
larger or smaller than lod- 

Now we evaluate the correction due to g^d ~ 9 ad- ^ n appendix C we show that an adiabatic solution of the motion 
equation of g R A is possible, in the sense that the functional dependence on to is the same as the equilibrium one, but 
the gap parameter changes slowly with time (see Eq. l|C2|l ): 



with 



C A) - +(-) 



M = 



M 



^(u,±tO+) 2 -\A(t)\< 



to A(t) 
-A(t)* -to 



(28) 



(29) 



This functional form for the R/A functions is obtained if the e — h symmetry is maintained and if one neglects the 
diffusion in space-time which will be mainly important at intermediate times^. 

This adiabatic approximation in the advanced and retarded Green's functions allows us to write the Keldysh propa- 
gator in the form: 



9 K = 9ad qmV2p 2{lo 2 - \A\ 2 ) ?/(V2co c ) 



To calculate f K we resort to the analogous of ea. l|Cl|l which is valid for g K : f K = g K A/uo. Hence we have: 

f K = fad - ?(*)4V2 - p \2(co 2 - |A| 2 ) V(V2wc)l for to » A . 



(30) 



(31) 



Now we insert Eq. (fTTTf> into Eq. (fT7)l and consider the linear correction to the gap of the irradiated contact according 
to : 



A(t) 



SA(t) 



Here SA(t) is the correction to the unperturbed gap parameter A due to the radiation. In the limit of zero temperature, 
up to first order in q(t) and A/u>d, SA(t) is given by: 



SA(t) 



ln(2co D /A) - 2 + 0(A 2 /uj 2 d ) JA 



dto 

— P 

to 



2{to 2 



\A\ 2 )?/{V2u; c 



(32) 



It is interesting to note that the correction arising from the adiabatic dynamics has the role of renormalizing the 
coupling q{t) via the prefactor — [/ji(2o>d/ A) — 2] _1 . This prefactor is negative because Iti(2lobI A) > 2. Therefore 
Eq. (|32|l shows that the gap of the irradiated superconductor is decreased due to the non equilibrium distribution of 
the qp excitations. 
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FIG. 2: Variation of the gap immediately after the pulse vs the inverse of duration of the pulse r c 



uj c in units of A. Our 



approximations break down for very low ratios w c /A (long pulses) 
this situation is unrealistic with the available optical devices. 



On the other hand, large oj c / A represent very short pulses: 



Eq. l|.{'2[l has the same structure as Eq. (14) of refill. There is a striking difference however. The inverse square 
root singularity at the gap threshold, which shows up in Eq.(14) of refill does not appear in Eq. 1|32|) . 

The inverse square root singularity originates from the unperturbed density of states at the excitation threshold 
and is usually present in non-stationary superconductivity^. It is responsible for retardation and oscillating tails. In 
our case, the gap threshold plays a little role, because we do not have extensive pair-breaking and q.p. generation at 
energies ~ A. Hence, just a tail ~ l/u> survives in the integrand. 

In Fig. (01 the variation of the gap immediately after the pulse (i ~ 0), is plotted vs the inverse of the pulse duration 
cu c in units of A, for different values of ojd- Our approximations are not valid when the pulse becomes too long (very 
low values of w c /A). For longer pulses the integrand in Eq. (|32fl has a narrow peak lined up at the gap threshold. 
In this case the inverse square root singularity in the density of states at the gap threshold is important and the 
adiabatic approximation Ea. 1)30(1 breaks down. 

For shorter pulses the peak becomes broader and is centered at larger w's. If the integration range is small, the 
result is quite sensitive to the location of the peak (see full line in Fig. (J2J): most remarkably, a minimum appears in 
the curve when the pulse is rather long (lo c /ujd < 5). By contrast, the gap correction is rather flat w.r.to changes of 
lu c when aj£>/A is larger (broken and dotted line in FigEJ. 



B. Correction to the Josephson current 



In this subsection we derive the correction to the Josephson current arising from the two terms of the anomalous 
propagator f(R,t,t') given by Eq. H31[). Ea.s H3UI31[) show that the non equilibrium Keldysh Green functions of the 



irradiated superconductor can be separated into two terms. The first one is what we call the 'adiabatic' contribution, 
while the second one is strongly dependent on the non equilibrium qp's distribution function and is first order in q(t). 
Within the linear response theory in the tunneling matrix element |T G | 2 , the pair current at zero voltage is: 

/oc 
dt , e 2 ie V(t-t') 
-oc 

(/ >f i(0,M') / 2 >,t',i) + f R l(0,t,t') f<(0,t',t)) , (33) 

where \T a \ 2 is assumed to be independent of energy, for simplicity. The current of Eq. (|33J) is evaluated at the junction 
site, defined by R = and the irradiated superconducting layer SI is labeled by 1 here, while the superconductor 
unexposed to radiation S2 is labeled by 2. The perturbed Josephson current has an adiabatic term J ad {t) obtained 
by inserting the first term of Eq. (|31J) into Eq. I|33|l . plus a correction 8J(t) arising from the second term of Eq. Ij31|l . 
Using the definitions f A = / > + / < and f R — f A = f y — f < and expanding to lowest order in q(t) the adiabatic 
critical current is: 

r d = J f{t) sinM, J fit) = ^ (l + ^) , (34) 

where Rn is the normal resistence, A is the unperturbed gap parameter of both contacts (assuming Ai = Ae lip and 
A2 = A in the absence of laser perturbation), and 5A(t)/ A is given by Eq. I|32[) . Denoting by Sf K the second term 
of Eq. I|31(l the correction SJ(t) is: 

6J(t) = e\T D \ 2 ^ duj ({8f K {Lo + ieV))\f£{uo - ieV) + (f R (u + ieV))\{8 f K \u - ieV))i) . (35) 

J — OC 

At zero temperature and V = 0, this gives: 

§J (t) = 4jL r dw^^-p\2(u; 2 -\A\ 2 )i/(V2cJe) 
eRw ' 



oo 



v / (w + i0+) 2 - |A| 2 ^(co-iO+y - |A| 2 



(36) 



which is zero for parity. This conclusion holds because we assume that no charge imbalance occurs. If V =/= the 
contributions to the integral evaluated in the complex plane are finite. 

8J(t) is a cosip-\ike correction. In the unperturbed Josephson effect a 'cosy' term only arises when V > 2Ai By 
contrast, our calculation shows that a costp term can arise in the Josephson current with a small nonzero voltage in 
the presence of an ultrafast laser pulse. 



IV. CLASSICAL DYNAMICS OF THE IRRADIATED JUNCTION 



In this section we integrate the classical equation of motion of the irradiated junction numerically. Here we discuss 
the possibility that the laser pulse induces switches of the junction from the zero voltage state, to the resistive 
state. The characteristics of the Josephson junction for a finite voltage, is obtained within the RCSJ (resistively and 
capacitively shunted junction) modeli. The phase of the superconductor S2 is taken as the reference phase. 

In the absence of the pulse, the junction is biased by a current constant in time As discussed in the previous 
Section, the pulse activates the superconductor SI by varying its gap dynamically in time. 

Consequently, a voltage V arises at the junction, related to the dynamics of the phase difference ip(t). The latter 
solves the differential equation: 

^ + %V+^TT i sin^)= 7 , (37) 

where 7 = Jb/J° an d J° = (whA) / (2cRn) is Josephson critical current of the unperturbed junction. The time- 
dependent driving term is deduced from Eq. J33J. We assume u> p jq > t^ 1 (r i ^ 1 (iV&) ~ 7GHz), where u> p jq is 
the plasma frequency at zero bias. This condition is satisfied for high quality Nb junctions, where lo p jo is in the 
range AQGHz ~ 120GHd^, but it holds also if we take into account the dependence of the plasma frequency on 7: 



Lo p j = uj p jq(1 — At 7 = 0.98 the term (1 — 7 2 )( 1 / 4 ) = 0.44: this still gives a large plasma frequency for the 

given range. In any case the plasma frequency changes marginally when the energy is degradated into heat if q is 
small. Under this conditions the relaxation process occurs long after the switch to resistive state. 

In Eq. I|37[) . Qo = u p jqR{lp)C is the quality factor, where R(ip) is the junction intrinsic resistance, which is in 
general a non-linear function of the phase. The dissipative Qq 1 (£ term includes thermal incoherent pair breaking 
effects at equilibrium. In the simulation we use both a constant junction resistance R and a patchwork model given 
by^: 

Q-\<p) = Q-l^£ I A ) (38) 

with N = 16 and Qq 1 — 0.636, which corresponds to a normal resistance above the gap Rn = (ui p jifo)/(cJ^). In 
general we ignore the direct dependence of R(<p) on the phase. By the way, Q should also depend on the energy which 
is released by each single laser pulse due to the incoherent pair breaking process. However, under the hypothesis 
that this energy is very low we assume that the quality factor, due to the optically induced normal resistance of the 
sample, is constant within the considered energy range. 

Actually, in the presence of the pulse, also the current contribution of Eq. I(35|) should be plugged into the l.h.s. of 
Eq. 137(1 . This current term depends on the voltage V = ip. However, in view of the fact that in this work we are 
only concerned with the switching of the junction out of the zero voltage state, we do not derive the full dynamics of 
the phase self-consistently. 

In Fig. @ we show the voltage just after the pulse for different values of the released energy. The time evolution 
of the voltage is sketched for some successfully induced switches. The junction starts in the zero voltage state. At 
Upjt = + it is irradiated by the laser pulse. There are few oscillations at frequency lu p j before the switching occurs, 
followed by an overall increase of the oscillating voltage. The larger is q the faster is the switch. If no switch is induced 
the junction remains in the zero voltage state: the phase and the voltage are weakly perturbed by the radiation and 
show decaying plasma oscillations around their equilibrium values. 

In Fig.(@J we show the approach to the gap voltage for different Q values and two different conductance models. 
Except for the asymptotic trend, the curves for Q = 10 (B) and Q = 100(A) show a similar behavior. The non-linear 
conductance gives rise to a more pronounced shoulder in the curve after the first increase of the voltage. The first 
phase oscillation at frequency oj p j are largely independent of the dissipation model used. 

The switching of the junction out of the zero voltage state depends on the bias current Jb, on the released energy 
per pulse g, and on the pulse duration. In Fig.© we sketch the switching front in the parameter space a)(w c /A, q) at 
fixed 7 and b) (7, q) at fixed uj c /A, for Q — 100 and wd/A = 10. For each point (q, w c /A, 7) we calculate J,? d (i) from 
Ea.s l|34l I32f> . Next we plug the result into the equation of motion Eq. (|37|l . Numerical simulation of the dynamics 
shows whether the junction is stable in the zero voltage state, or it switches to a running state. The points of the 
curve mark the frontier between the two behaviors. The full curve is just a guide for the eye. The non-monotonicity of 
<5A(t)/A with the pulse duration, appearing in Fig- (J2J forces a similar behavior in Fig.|J5Ji). This means that the pulse 
duration can be appropriately chosen, in order to optimize the junction switching with the laser field. Indeed, if uj c /A 
is of the order of 5 -j- 15 a very small released energy is required for the switching of the junction, because the order 
parameter is much depressed by the laser pulse in that range. Out of this range the shorter is the pulse the larger is 
the energy required. By contrast a slightly larger released energy is also required for longer pulses (u) c /A < 5). This 
is because longer pulses imply a more extended change in the qp distribution up to higher energies. As a consequence 
the maximum of the function p of Eq. |22contributes less to the correction of the gap parameter. Nevertheless caution 
should be used in considering our results for longer pulses because of the neglected relaxation effects. 



V. CONCLUSION 



The effect of an ultrafast laser pulse on the superconducting coherence at a Josephson Junction allows studying an 
unexplored regime in non equilibrium superconductivity. 

Non-equilibrium in superconductivity is usually addressed in the context of one of the possible applications of 
Josephson junctions, that is radiation/particle detectors. Highly energetic radiation produces pair breaking and 
quasiparticles which, in turn, excite a large number of them, in a cascade process. Usually the setup is optimized such 
that the qp's can be collected and contribute to the current across the junction with a sharp signal. Losses are due 
to degradation of the released energy into heat during the relaxation process. To achieve optimum performance, the 
Josephson current is usually suppressed by applying a magnetic field. This picture has been discussed quantitatively 
by examining the quasiclassical kinetic equation for the non equilibrium qp distribution function^. 




FIG. 3: Voltage behavior in time for different energy released on the sample and the quality factor Q = 100. The voltage is 
normalized to Vo = u p jifo/c- 

In this work we have concentrated on a quite different time scale: the one fixed by the duration of an ultrafast laser 
pulse. While the relaxation process mentioned above takes place on a time scale of 0.1 -j- 100 ns, we have considered 
a laser perturbation lasting at most hundreds of femtoseconds. This type of tool can be quite valuable for future 
applications, because fast pulses in flux and gate voltages are extremely important when processing information in 
superconducting quantum computing devices (qubits)2i. Indeed, finite rise and fall times of pulses may result in a 
significant error in dynamical computation schemes^. The carrier frequency of the laser is ~ lOOTHz and the optical 
radiation is expected to produce many e — h pair excitations as it would occur in a normal metal. In our case, qp's 
do not have enough time to relax down to the typical phonon frequencies (~ ljd) and to heat the sample before the 
stimulated switching occurs. We do not wish to collect qp's either, what requires a suitable geometry of the junction. 

Instead, we have addressed the question how the critical current for Josephson conduction J c can be coherently 
affected by a laser induced small perturbation with an unrelaxed non-equilibrium distribution of qp's, that is before 
the dissipative response sets in. 

Using the quasiclassical approach to non-equilibrium Keldysh Green's functions, we have shown that, if the temper- 
ature is very low, the order parameter of the irradiated superconductor can respond adiabatically to a weak perturbing 




FIG. 4: (Color online) Voltage behavior in time with different quality factors Q in the linear conductance model (A,B) and the 
non-linear conductance model as given by Eq. I|M8^ (C). The voltage is normalized to Vo = ui p jipo/c- 



signal. A non equilibrium distribution of qp's is generated and consequently J c is temporarily reduced (see Eq. 1|34|) 
and FigQ). This reduction can drive switches of the junction out of the zero voltage state. In our approach the 
retardation effects which arise from the frequency dependence of the e — ph coupling a 2 F(u>) and from the actual 
features of e — ph relaxation processes have been neglected. They came into play on a time scale longer then the 
duration of the laser pulse, lu^ 1 ~ r c . Indeed in the equation of motion for the Keldysh Green's function reported in 
appendix B and C the role of the frequency dependent self-energy terms has not been discussed. 

The parameter which is related to the energy released by the radiation and describes the strength of the perturbation 
is q. In our case, the switches can be induced by pulses with q ~ 0.05, with relatively low values of uj c /A ~ 10, by 
polarizing the junction very close to the critical current. 

In experiments on laser induced non-equilibrium effects in superconductors- 2 or Josephson junctions^, the released 
energy is of few fi J, which corresponds to values of q between 0.25 and 0.67 for the given laser spot dimension. In 
our case for q = 0.05 the switching occurs at 98% of the critical current. Therefore, a coherent effect of the laser 
on the superconducting condensate is sufficiently large to be observed in sensitive experiments monitoring the escape 
rat o 1 ^ 11 . These experiments can appreciate very small variations of the critical current, if temperature is kept low 
and the released energy is sufficiently small, so that sizeable heating effects do not occur. 

We have also found a 'cosy' contribution to the Josephson conduction due to the presence of the excited qp's 
(see Eq. (j2HJ)- This term, which will be examined in detail elsewhere, vanishes at zero V,T, provided excitations 
do not generate charge imbalance. A similar term, proportional to the voltage V, can be derived also in the BCS 
theory of Josephson conduction^ but it is identically zero as long as V < 2A/e, because of the absence of qp's at zero 
temperature. This is not the case here, due to the presence of a non-equilibrium qp distribution. 

Our derivation of Ea.s (|32l34l36J) assumes that no charge imbalance is created by the perturbation. This is because 
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FIG. 5: Switching front in the parameter space a)(w c /A, q) at fixed current bias 7 = Jb/J° an d b) (7, q) at fixed cj c /A, for 
Q — 100 and wd/A = 10. The full curves have been added as a guide for the eye: they mark the border between the zero 
voltage (Josephson) state and the voltage state, q is the coupling strength due to the laser pulse. 



the radiation excites e — h pairs and the pulse duration is short enough, so that pair breaking is very limited. The 
absence of charge imbalance is a crucial approximation in our solution scheme. This assumption allows us to keep 
the unperturbed functional form of the quasiclassical Green's functions and to insert a time dependent gap parameter 
A(£) in their expressions. A(i) follows the perturbation adiabatically and is determined by the non equilibrium e — h 
pair distribution produced by the pulse. Charge imbalance corrections should be reconsidered, but they are usually 
expected to have a minor effect on the dynamics of the junction. 

To complete the picture, we have simulated the classical dynamics of the junction switching to the resistive state. 
This picture is only valid over few periods 2tt /co p j on time duration less than the electron- phonon relaxation time. 

Precursor oscillations can be seen in Fig. J2J at the Josephson plasma frequency u> p j. Most remarkably the duration 
of the pulse can be optimized in order to induce controlled coherent switching at the minimum possible released energy 
£. 
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APPENDIX A: QUASICLASSICAL TIME DEPENDENT GREEN'S FUNCTION APPROACH 



The quasiclassical Green's function solve the Eilenberger form of the Gorkov equations in commutator formal: 



G^-E-A ,G 



= 



(Al) 



The matrix Green's function G in the Keldysh space is: 



G = 



G R G K 
6 G A 



where, in turn, G R ,G A ,G K are the retarded, advanced and Keldysh Green functions in the Nambu space: 



q(A,R,K) 



g(A,R,K) j-(A,R,K) 
_f(A,R,K)i _ g {A,R,K)\ 



Here / is the anomalous propagator and its Keldysh component defines the gap: 

z/(0)A r+ °° 



A 



duj < f K 



The average <>v F denotes angular averaging over the Fermi surface. The gap matrix A is defined as: 



A = 



A 
A 



A = 





A* 



-A 




(A2) 



(A3) 



(A4) 



(A5) 



The self-energy S includes elastic and inelastic scattering with impurities and gives rise to relaxation processes. The 
commutator is evaluated w.r. to the (£> operation, which implies integration over the intermediate variables according 
to: 



Gn 1 ® Go 



d2 Go 1 (l,2)G(2,l') , 



where 1 = (ri,ii). The differential operator G (1,2) is 



<V(M) 



i~cM(l) +( e 0(l)- M )/ 
c 



6(1 2) 



Here 



o t 
&i 



I 



I 

o i 



(A6) 



(A7) 



(A8) 



where &i(i = 1,2, 3) are the usual Pauli matrices and / the 2x2 unit matrix 
The vector potential ^4(1) describes the laser radiation field of frequency O: 



=pz(p-r*i — Ot) 



(A9) 



where a± can be slowly varying 'envelope' functions of space, on the light spot size R and on time, on the scale 
n _1 . These are the reference space and time scales in the following. In the frame of the quasiclassical approximation, 
the original Green's functions (2(1,2) = G^^,^ ,t') are assumed to be slowly varying function of the coordinate 
R = (r + r*)/2 while they oscillate fast as functions of r — f on the scale of the Fermi wavelength Xp- It is customary 
to rewrite also the time dependence in terms of the new variables t = (t + 1')/2 and t — t' and to Fourier transform 
w.r.to f — f and t — t', thus obtaining G(p,u),R,t). 
The motion equation for the Keldysh component of ea. ljAlj) reads: 



G ( 7 1 -Ke(S)-A,G 



K 



Z K 7 3le(G) 



-\± k .AX --(f, G K \ 



(A10) 



where we have defined a quantity proportional to the density of states A — i(G R — G A ) (not to be confused with 
the vector potential), and written down the imaginary and the real part of the self-energy, T — i(T, R — T, A ) and 
JieE = J; (S^+S^), respectively, as well as the real part of the retarded/advanced Green's function SfteG = ^(G R +G A ) 
({, } denotes the anticommutator ). The next step is the gradient expansion of the £g) product (we drop the overline 
on t in the following): 



C ® B = exp l/2(« - ) exp 1/2(5^^ - d^)CB 



(All) 



(here dp stands for the gradient w.r.to the impulse operating on C ), followed by the averaging of the result for \p\ 
close to pf, that is over the energies p 2 /2m — ^ = £ while the direction of p, j3, is untouched: 

g(p, Lu,R 1 t)= % - J dt G(p, w, R, t) . (A12) 

This is justified, because is much shorter both of the superconducting correlation length and of the spatial range 
of the laser spot (e — h symmetry is assumed). Eventually g depends on p,oj,R,t. The average of ea. ljAlOjl over 
all directions in the Fermi surface, can be done if no external bias is applied and anisotropies of the diffusion and 
relaxation process are not expected. In the presence of radiation with optical frequency it is customary to average 
out the fast oscillating components with frequency fi 28 . Following cq.JSJ), we expand the Green's functions similarly: 

g(uj, R,t)=g(Lu,R,t)+J2 X ± (w, R, t)e ±mt ■ (A13) 

± 

Here g(uj,R, t) is assumed to be the slowly varying part on the scale of the pulse duration r c , while are fast varying 
ones. All these functions are slowly varying functions of space as well, on the light spot size scale. A decomposition of 
Eq. IjAlOjl into harmonics arises. We are interested in the zero order harmonic equation, which shows a slow dynamics 
that can be followed coherently by the irradiated superconductor. Leaving the self-energy terms in Eq. l)A10jl for the 
moment aside and dropping the superscript ( K \ we obtain 

[na+, x~] + [r 3 a-, x + ] + [(wfa - A),£] - i{f 3 , d t g) + i{d t A, d u g) - ^^9 t A 2 d u g = ... , (A14) 

where the ellipsis refers to the missing self energy terms. The first order harmonic equations are: 

±i{T 3)X ± } + ^^[a±,g]=Q ) (A15) 

they show that the first two terms in Eq. (|A14|I are 0(f2 -1 ) smaller than the others and can be neglected to lowest 
order. Hence the effective equations for the Green functions are 

[{wh ~ A), g] - \{hid t g} + l{d t A,d u g} - ^0^8^=... , (A16) 

The last three terms in Eq. I|A16(1 include the time dependent non equilibrium dynamics that is absent in the case of a 
time independent approach. Retarded, advanced and Keldysh Green's function, they all satisfy analogous equations. 



APPENDIX B: KINETIC EQUATION FOR n(w,r,t) 

One can linearize Eq. (|A16J) for Keldysh Green's function by posing g K — g R h — hg A . This yields the kinetic 
equation for the distribution function n(u,R,t)^. We neglect any variation in space and concentrate on the t 
dependence here. The distribution matrix h is defined starting from the and nx functions according to: 



n L l + n T a- 3 



(Bl) 



or 



h = 



n L {E)+n T (E) \ = fl-2n(E) 

n L (E) - n T (E) J ~ I 2n{-E) - 1 



(B2) 



The second equality defines the relation with the qp distribution function n(E). We always assume e — h symmetry, 
so that n(E) + n(—E) = 1 and tit = 0. In the equilibrium case one has n {E) — e E/r +1 , so that: 



l L(T) 



E E 
ta nh( — ) + (-)tanh(-] 



(B3) 



Substituting Eq. I|19|) m Eq- (|A16(1 . we get the kinetic equation for the distribution matrix h.ln particular, in case 
there is no charge imbalance, the equation for the longitudinal component til is~: 



d t n L Tr((g R h - T 3 g A )) + d u n L Tr(d t A(g R - g A )) - —d t A 2 Tr (d u n L (g R - g A )) = -4I[n L (u)] 



(B4) 



where I[til{uj)] is a collision integral which regulates the relaxation of the qp distribution toward equilibrium. Eq. 
(|B4|I is fully discussed in refw 



APPENDIX C: EQUATION OF MOTION OF G 



We now write down the equations Eq. (|A16|) explicitly for the retarded Green's functions. We label the matrix 
components by — 1,2) and drop the superscript R everywhere. The matrix elements of Eq. ljA14j) in the 

Nambu space become : 



(1, 1) -> d t g - A/t + A*/ + d t Ad u p + d t A*d^f + -^-d t A 2 d u g = ... 

Zmc z 

(2, 2) 9 t5 t + A/t - A*/ + ft Ad„p + d t A*d w f - -^-ft^cW = ... 

2toc z 

(1, 2) 2w/ - A.gt - gA + ftAcW _ 9w59tA + JL_Q tA *Q^ = ... 

(2, 1) +2c^/t - A* 5 - g^A* + d t A*d u g - d^d t A* - ^9 t A 2 d u p = 



The equilibrium result suggests that 



S-*&. (CD 



solves Eq.(l, 2) except for terms oc dtA 2 which describe the relaxation at later times. Let us assume that this relation 
holds also in the non-equilibrium case. Then the formal solution, follows adiabatically the t— dependence of the gap 
parameter A by keeping an equilibrium-like shape: 

9ad ~ V^-IAPW' " v^ 2 -|A| 2 W [ ' 

This approximate solution is quite appealing, because it satisfies the equilibrium condition for A at t — > oo. However 
it neglects diffusion in space-time which will be mainly important at later times w.r. to the pulse duration. 



1 A. Barone and G. Paterno, Physics and applications of Josephson effect, Wiley New York, 1982. 

2 L. R. Testardi, Phys. Rev. B 4, 2189, 1971. 

3 W. H. Parker and W. D. Williams, Phys. Rev. Lett. 29, 924, 1972. 

4 C. S. Owen and D. J. Scalapino, Phys. Rev. Lett. 28, 1559, 1972. 

5 R. A. Vardanyan and B. I. Ivlev, Sov. Phys. JEPT 38, 1156, 1974. 

6 A.V. Sergeev, M.Yu. Reizer, J. Mod. Phys. B 6 635 (1996), M.Yu. Reizer, Phys. Rev. B 39 1602 (1989), A.V. Sergeev, 
M.Yu. Reizer, Zh. Eks. Teor. Fiz. 90 1096 (1986) [Sov. Phys. JETP 63 616 (1986)] 

7 M. Lindgren, M. Currie, C. Williams, T.Y. Hsiang, P.M. Fauchet, R. Sobolewski, S.H. Moffat, R.A. Hughes, J.S. Preston, 
and FA. Hegmann Appl. Phys. Lett. 74 853 (1999) 

s A.D. Semenov, R.S. Nebosis, Yu. Gousev, MA. Heusinger, K.F. Renk, Phys. Rev. B 52 581 (1995). 

9 S.B.Kaplan, C.C.Chi, D.N.Langenberg, J.J.Chuang, S.Jafarey, D.J. Scalapino, Phys.Reb.B 14, 4854, 1976. 

10 M. H. Devoret, J. M. Martinis, D. Esteve and J. Clarke, Phys. Rev. Lett. 53, 1260, 1984. 

11 M. H. Devoret, J. M. Martinis and J. Clarke, Phys. Rev. Lett. 55, 1908, 1985. 

12 A.Moopenn, E.R. Arambula, M.J. Lewis, H.W. Chan, IEEE Trans. Appl. Supercond. 3, 2698, 1993, 

13 G. P. Pepe, G. Peluso, M.Valentino, A. Barone, L. Parlato, E. Esposito, C. Granata, M. Russo, C. De Leo and G. Rotoli 
Appl. Phys. Lett. 79, 2770, 2001. 

14 Superconductive Particle Detectors, edited by A. Barone (World, Singapore, 1988); A. Barone, Nucl.Phys. B 44 645 (1995). 

15 Yu.N.Ovchinnikov and V.Z.Kresin Eur.Phys.J.B 32, 297, (2003), Yu.N.Ovchinnikov and V.Z.Kresin, Phys.Rev. B. 
58,12416,(1998). 

16 A.I.Larkin and Yu.N.Ovchinnikov Zh.Eksp.Teor.Fiz 55,2262 (1968) [Sov.Phys. JETP 28, 1200, 1969]; J. Rammer, H.Smith, 
Rev. Mod. Phys., 58, 323,(1986). 

17 W.Belzig, F.K.Wilhelm, C.Bruder, G.Schoen, E.D.Zaikin Superlattices and Microstructures 25, 1251 (1999). 

18 M.H.S.Amin, Phys.Rev.B 68 , 054505, (2003). 

19 preliminary report on some of this material has appeared on P.Lucignano, F.J.W.Hekking and A.Tagliacozzo, in "Quantum 
computing and Quantum bits in mesoscopic systems", A.J.Leggett, B.Ruggero, P.Sivestrini ed.s, Kluwer Academic, Plenum 
Publishers N.Y. (2004). 

20 G.Grimvall 'The electron phonon interaction in metals', North Holland Publ. Co. Amsterdam 1981 

21 J.J.Chang and D.J. Scalapino, Phys.Rev.B 15 , 2651, (1977). 



L.P.Gorkov and G.M.Eliashb erg, Soviet Physics,JEPT27,328 (1968) 

A.A.Abrikosov, L.P. Gor'kov, I.E.Dzyaloshinski, Methods of quantum field theory in statistical physics; 

A. Wallraff, A. Franz, A. V. Ustinov, V. V. Kurin, I. A. Shereshevsky and N. K. Vdovicheva, PhysicaB 284-288, 575, 2000. 
K. K. Likharev, Dynamics of Josephson junctions and circuits, Gordon and Breach Science Publishers, 1986. 
Yu. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. Phys. 73, 357 (2001) 
M.S. Choi, R. Fazio, J. Siewert, and C. Bruder, Europhys. Lett. 53, 251 (2001) 
L.D. Landau, E.M.Lifsits, Mechanics, Pergamon Press (1960). 



